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1 Abstract 



We compute the value of a variance swap when the underlying is modeled as a Markov process 
time changed by a Levy subordinator. In this framework, the underlying may exhibit jumps with a 
state-dependent Levy measure, local stochastic volatility and have a local stochastic default intensity. 
Moreover, the Levy subordinator that drives the underlying can be obtained directly by observing Eu- 
ropean call/put prices. To illustrate our general framework, we provide an explicit formula for the value 
of a variance swap when the diffusion is modeled as (i) a Levy subordinated geometric Brownian motio n 
with default and (ii) a Levy sub ordinated Jump-to-defau l t CEV process fsee lCarr and Linetskvl (|2006h ). 
Our results extend the results of iMendoza-Arriaga et al.l (|201Cf ). by allowing for joint valuation of credit 
and equity derivatives as well as variance swaps. 

1 Introduction 

A variance swap (VS) is a forward contract written on the realized variance of an underlying S = {St, t > 0}. 
As is typical in derivatives literature, we define the realized variance over the interval [0,f] as [log S]t, the 
continuously sampled quadratic variation of the log S. Thus, at maturity the VS has a payoff (to the long 
side) of 

[log S] t -K var . (1.1) 

The variance swap rate K var = E ([log S]t) is determined at inception so that initial value of the VS is zero. 
Here, E is the expectation under the risk-neutral pricing measure P. 

There are a number of reasons for which one may wish to enter into a VS agreement. First, a trader 
who delta-hedges a short p osition in a European optio n can limit his exposure to stochastic volatility risk 
by trading a VS (see. e.g. Carr and Schoutens ( 2008| )). Second, a drop in the level of an underlying S is 



often accompanied by an increase in the volatility of S (the leverage effect). Thus, a long pos ition i n a V S 



serves as protection against a market crash. Such is the demand for VSs that, according to I Jung (|2006h . 
the daily trading volume in equity index VSs reached 45 million USD vega notional in 2006 (vega measures 
the change in an option's price caused by changes in volatility). On an a nnual basis, th i s corr esponds to 



payments of more than 1 billion USD per percentage point of volatility (see Carr and Led ( 2009h ) 



At the individual stock level there is yet another source for the leverage effect that is related to the 
default risk associated with the underlying firm. This credit-related l everage effect explain s the interaction 
between market risk (return variance) and credit risk (default arrival). ICarr and Wu ( 2009t ) recently studied 



this interaction in pricing stock options and credit default swaps. In general, when the probability of default 
of a firm increases, its stock price tends to lose value while the option-implied volatility increases. In such a 
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situation, it is common for investors to look for credit protection by taking long positions in deep out-of-the- 
money puts. However, there is a side-effect in this hedging strategy. Though the investor's objective is to 
minimize credit exposure, this strategy will increase his exposure to volatility due to the credit-equity leverage 
effect. Hence, in the pricing of variance swaps written on an individual stock it is especially important to 
take into account defaul t risk (and t he cr edit-related leverage effect). 

In his seminal paper, Neuberger ( 1990f ) showed that, when the underlying S is modeled by a process with 
continuous sample paths, dSt = cr t StdWt, (for simplicity, assume the risk-free rate of intere st is zero) the fair 
value of K var is given by a European-style log contract K var = E[— 2 log(S t / Sq)]. Later. ICarr and Madan 
(1998) showed that the log contract (or in fact, any European-style derivative with a twice differentiable 
payoff) can be synthesized from a linear combination of calls and puts. Q Thus, when i-expiry calls and/or 
puts are available at every K € (0, oo), the fair value of K var is uniquely determined by the implied volatility 
smile. Remarkably, this result is independent of any assumption about the volatility process a = {o- t , t > 0}. 
However, this result does rely on the continuity of sample paths of S. The events of the recent (and ongoing) 
financial crisis underscore the need to include jumps in the underlying price process S. The question naturally 
arises then: what is the fair value of K var when the price process S is allowed to have discontinuous sample 
paths? 

One possible answer to this question is given bv lCarr et al. ( 2012 ). who show that, when log S is modeled 
as a Levy process time-changed by an absolutely continuous time-change, the fair value of K var is equal to 
a multiplier Qx times the value of a log contract K var = Qa'E [— log (St/ So)]- Interestingly, the multiplier 
Qx does not depend on the time-change process. While the value of K var does depend on the time-change 
through the log contract, since the log contract can be synthesized by a linear combination of calls, when 
one has fu l l know ledge of the volatility smile, the value of K var can be determined in the framework of 
Carr et al.l ([20121 ) without any knowledge of the time-change process. 

This last point cannot be emphasized enough. A parametric model for the time-change process would 
leave one open to model risk, as any misspecification of the time-change parameters would (in general) result 
in errone ous values of K va r . By using knowledge of the volatility smile to construct the value of the log 
contract. ICarr et al.l (|2012f) circumvent the need to parametrically model the time-change process. Thus, the 
risk of model misspecification is greatly reduced (though, some model risk still exists, since the multiplier Qx 
depends on the choice of a specific Levy process). An alternative and quite interesting a pproach to the joint 
pricin g of volatility derivatives and index options in the presence of jumps is provided in lCont and Kokholm 
1)201 lh . 

In this paper, we consider the class of Levy subordinate d diffusion process es described in lMendoza-Arriaga et al 
( 2010 ). This class of models, like the class considered by ICarr et al. ( 2012 ). allows for the underlying S to 
experience both jumps an d stochasti c volat ility. However, there are a few important differences between 
these two frameworks. In ICarr et al.l ( 20121 ). the background process is modeled as a Levy process, which 
naturally includes the possibility of jumps, but does not include stochastic volatility. Stochastic volatility is 
a dded by time-changing the bac kground process with a continuous increasing stochastic clock. In contrast, 
in lMendoza-Arriaga et al". (2010). the background process is modeled as a diffusion, which may include local 
stochastic volatility, but does not include jumps. Jumps are added by time-changing the ba ckground pro- 
cess with a Levy subordinator. Additionally, the frame work of Mendoza-Arriaga et al. ( 20101) allows for the 
possibility of a default event, whereas the framework of Carr et al.l (|2012l) does not. While default may not 
be realistic for an index, it is certainly an important consideration for individual stocks as explained above@ 

Despite the differences between these frameworks, they share one desirable feature: when the background 
process is fixed and one has full knowle dge of the volatility smile, the fair value of K var is robust to misspec- 
ification of the time-change process. In lCarr et al.l (|2012l ). the effect that the time-change has on the value 
of K var is felt through the log co ntract, which is constructed directly from European calls. In the framework 
of lMendoza-Arriaga et al.1 ()2010h . we will show that the Levy subordinator can be inferred directly from the 
volatility smile. Once the subordinator is obtained, it can be used to compute the fair value of K var (among 



1 In fact, the replication result of lCarr and Madanl l|l99 8f) is independent of the continuity assumption of the price process. 
2 A default event would cause [log S]t to blow up. As such, we must amend our definition of a VS to account for this 
possibility. We will do this in section [4] 
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other things). 

The rest of this paper proceeds as follows. In section[2]we describe the class of Levy subordinated diffusion 
processes in detail. In section [3] we introduce some mathematical tools, which we shall need to compute the 
price of a VS. In section 2] we modify the payoff of a VS to account for the possibility that the underlying 
defaults (i.e., jumps to zero). We then derive a general expression for the value of the modified VS in the 
Levy subordinated diffusion setting. In section [5] we present some important results concerning generalized 
eigenfunction expansions. These results will be needed for the diffusion-specific VS computations provided 
in section[B] In section[7]we show that, by observing European call and put prices written on the underlying, 
one can uniquely determine the drift and Levy measure of the subordinator that drives the underlying price 
process. Finally, in section [51 we provide some concluding remarks and discuss directions for future research. 



2 Model 

We assume frictionless markets, no arbitrage, and take an equivalent martingale measure (EMM) P chosen 
by the market on a complete filtered probability space (f2, 5F, {3^, t > 0}, P) as given. All stochastic processes 
defined below live on this probability space, and all expectations are with respect to P unless stated otherwise. 
We also assume that all the (natural) nitrations defined throughout this section are augmented to include 
the P-null sets. In particular, this assumption holds for the subordinate filtration £Fj, which we describe at 
the end of this section. For simplicity, we assume zero interest rates. All of our results can be easily extended 
to includ e deterministic interest rates. 

As in Mendoza-Arriaga et al. ( 201dt ). we model the stock price dynamics under the risk- neutral pricing 



measure P as a stochastic process S — {S t ,t > 0} defined by 

S t = (l-D Tt )eP t X Tt , A =%<«}, X = x, T f = 0, C>0. (2.1) 

Here, the background process X = {X t ,t > 0} is a scalar Feller diffusion, T = {T t ,t > 0} is a Levy 
subordinator, p is a scaling factor, which is needed to ensure that the asset price S is a martingale, and £ is 
a positive random variable, which will be used to model a possible default event of the underlying S. Below, 
we describe each of the above-mentioned elements in detail. Additionally, we discuss how some important 
nitrations relate to a trader's observation of the market. 

Background Feller process X. We let X = {X t ,t > 0} be a time-homogeneous Markov diffusion 
process, starting from a positive value Xq — x > 0, which solves a stochastic differential equation (SDE) of 
the form 



where 



dX t = b(X t ) dt + a(X t ) dB t , (2.2) 



b(x) := \pt + k(x)]x and a(x) :— cr(x)x. (2.3) 



Here, o~(x) and [/j, + k(x)} are the state-dependent instantaneous volatility and drift rate, /i € M is a constant, 
and B = {B t ,t > 0} is a standard Brownian motion. We assume that o~(x) > and k{x) > are Lipschitz 
continuous on the interval [e, oo) for each e > (i.e., locally Lipschitz), and that o~{x) and k{x) remain 
bounded as i — > oo. We do not assume that o~(x) and k(x) remain bounded as x — > 0. Under these 
assumptions the process X does no t explo de to infinity (i.e., infinity is a natural boundary for the diffusion 
process; see Borodin and Salminenl ( 20021 ) p. 14 for boundary classification of diffusion processes). We also 



assume that zero is either an (unattainable) natural boundary or an entrance boundary. If zero is a natural 
boundary the state space is given by E = (0, oo). If zero is an entrance boundary, i.e., the process X can 
be started at x = then it quickly moves to interior of [0, oo) to never hit zero again. Throughout this 
document we assume that the process always starts from a positive value Xq — x > 0, and hence the state 
space is also defined as E = (0, oo). Under all our previous assumptions X is the unique strong solution to 
the SDE (|2.2[) . The transition function P^{x,dy) = F(X t € dy\X n = x) of the diffusion process X started 
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at x defines a Feller semigroup CP = {7^,t > 0} acting on the space C([0, oo]) of functions continuous on 
(0,oo) and such that the limits lim^^o f(x) and lim^-j.oo f(x) exist and are finite (Ethier and Kurtz (1986) 
p. 366) by 

?°J(x)=E x [f(X t )} = [ f(y)P?(x,dy). (2.4) 



The infinitesimal generator of T is a second-order differential operator of the form 

1 d 2 f df 
A°f(x) = -a 2 (x)^(x)+b(x)£(x) 

with the domain Dom(A°) = {/ £ C([0,oo]) n C 2 ((0, 00)), A°f £ C([0,oo])} if zero is an inaccessible 
boundary (natural or entrance). We also note that the semigroup leaves the space Cq((0, 00)) C C([0,oo]) 
of functions continuous on (0, 00) and having zero limits lim^^o f( x ) — and lim K _ i . 00 f(x) — invariant 
and is a Feller semigroup on it. Lastly, we denote by J x the natural filtration of the process X. 

The trigger event ( and the indicator process D. Let £ ~ Exp(l) be an exponential random 
variable, independent of X, and define the trigger event time Q as 



C = inf jt > : J k(X s )ds > £ j. 



That is, £ is the first jump time of a doubly stochastic Poisson process with jump intensity given by the 
killing rate k(x). Observe that the killing rate is added in the drift (|2.3p to compensate for the killing (jump- 
to-default). This compensation will be needed in order to ensure that the stock price S is a martingale. 
Moreover, since the process X cannot reach zero from the interior of the state space, then the underlying 
stock price process S cannot go to zero continuously, rather, it may only jump to zero from a strictly positive 
value. We denote by CP 1 = {CPj , t > 0} the Feynman-Kac semigroup associated with the killing rate k(X): 

V 1 J(x)=E 3: [e-ti k ^ du f(X t )] = / f(y)P?(x,dy), t > 0. (2.5) 



E 



V 1 is a (sub-Markovian) Feller semigroup on C([0,oo]) with the generator 

A 1 f(x)=A°f(x)-k(x)f(x) (2.6) 

with the domain Dom^ 1 ) C Dom(.A ). More precisely, Dom^ 1 ) = {/ £ C*([0,oo]) n C 2 ((0, 00)), .A 1 / £ 
C([0, 00])}. Since zero is an inaccessible boundary (natural or entrance) of the diffusion with killing at the 
rate k(x) it suffices to restrict the domain to Dom(yl 1 ) = {/ £ C n ( f0, oo]) D C 2 ((0, 00)), A l J £ C n([0,oo1)| 
whenever we start the process from the interior of E, i.e., x > fcf. iBorodin and Salminenl (|2002h pp.l5ff). 

Since the random variable C is not H x measurable, we introduce an indicator process D = {D t , t > 0} in 
order to keep track of the event {£ < t}. The indicator process D is defined by 

D t =%<*}, (1- A) =%>*}■ (2.7) 

Lastly, we denote by 5" D to the natural filtration of th e process D. 

The auxiliary bivariate process [X, D). As in iMendoza-Arriaga and Linetskvl ( 2012 ) we define the 



auxiliary bivariate process (X, D) = {(X t ,D t ), t > 0} with state space E = E x {0,1}. The process (X, D) 
is a Feller semimartingale in the enlarged filtration J X ' D = J* V 3 rD (3 ' D is the smallest filtration that 
contains J x and in which £ is a stopping time). Moreover, any function f(x,d) £ C([0,oo] x {0,1}) is 
represented as, 

f{x,d) = u(x)(l-d)+v(x)d=(u-v)(x)(l-d) + v(x), u,v £ C([0, 00]). (2.8) 

As we shall see below, the function u(x) = f(x,0) can be interpreted as a promised payoff function if the 
triggering event £ does not occur by time t > 0, while v(x) — f(x,l) can be understood as a recovery payoff 
function if the triggering event occurs prior time t > 0. 

The following theorem gives the Markovian characterization of the bivariate process (X, D). 
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Theorem 2.1. (i) The bi-variate process (X,D) is a Feller process with the Feller semigroup 7 = {CP f , i > 0} 
acting on f £ C([0, oo] x {0,1}) according to: 

y t f(x,d) = y° t v(x) + (1 - d)y\{u - v)(x), 

where u, v £ C([0,oo]), 7° is the Feller semigroup of the process X (|2.4[) and 7 1 is the corresponding 
Feynman-Kac semigroup (|2.5[) . 

(ii) The infinitesimal generator of the Feller semigroup 7 is given by 

Af(x, d) = A°f(x, d) + k(x)(f(x, 1) - f(x, d)) 

= A°v(x) + (1 - d)A 1 (u - v)(x), 

where A and A 1 are the generators of 7° and 7 1 , respectively. 

(Hi) If f(x,d) £ Dom(yi) (i.e., f is of the form (|2.8I) with u,v £ Dom(.A 1 ) y ) and (X,D) is the bi-variate 
process with Xq = x > and D$ = d £ {0, 1} ; then the process 



M{ := f(X t ,D t ) - f(x,d) - f Af(X s ,D s )ds 

Jo 



is a 3 ' -martingale. 



Proof. The proof can be found in iMendoza-Arriaga and Linetskvl (|2012j). 

Similarly, the Doob-Meyer decomposition of D is given in the following corollary. 
Corollary 2.2. The increasing process D has the Doob-Meyer decomposition 

D t = M t + A t 

with the predictable $ X ' D -compensator At and 1 X,D -martingale Mt 



□ 



A, 



f k(X u )(l - D u )du, M t = D t - [ k{X u ){ 
Jo Jo 



1 — D u )du. 



Proof. See secti o n 5.1 .4 in Bielecki and Rutkowski ( 20041 ) (also, Lemma 7.3.4.3 and Comment 7.3.4.4 in 
Jeanblanc etaH (|2009h ). □ 



The Levy subordinator T. A Levy subordinator T — {T t ,t > 0} is a Levy process with positive j umps 
and non- negative drift. For a standard reference on subordinators we refer the reader to lBertoinl (J200J). We 
require that T be independent of X and £ and satisfy To = 0. Every Levy subordinator has the following 
ltd- Levy decomposition 



dTt = 7 dt ■ 



sdN t (ds), 



(o, 



where 7 > is the drift of the subordinator and N is a Poisson random measure with the property that, 
for any Borel set A £ 23((0, 00)) we have E,[dNt(A)] — v{A)dt for some cr-finite measure v on 23((0,oo)). 
The measure u, which must satisfy J (s A 1) v{ds) < 00, is referred to as the Levy measure. The Laplace 
transform of a Levy subordinator T is given by 



E [e ATt ] = f 
Jo 



e Xs TT t (ds) 



X £ < X £ 



[1,00) 



e Xs v(ds) < 00 > =: J, 



(2.9) 
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where nt(ds) is the transition function of the subordinator T and 4>(X) is the Laplace exponent of T, which 
can be computed explicitly from the Levy-Kintchine formula 



0(A) = 7 A + / (1 - e~ As ) v{ds). (2.10) 



Note that the Laplace exponent <j>(X) is concave and increasing and satisfies 0(0) = fsee iBertoinl (|1998h . 



page 73). For any Borel set A e S((0, oo)) the process {N t (A),t > 0} is a Poisson process with intensity 
v(A) = J A v(ds). If the arrival rate of all jumps is finite a :— v((0,oo)) < oo then J Q sdN t (ds) is a 
compound Poisson process. In this case, the distribution of jumps is aT x v. If z/((0,oo)) = oo then the 
subordinator is said be an infinite activity subordinator. 

We define L = {L t , t > 0} the first passage process or right inverse process of T as L t := inf{s : T s > t}. 
Recall that T is assumed to be independent of S^' 15 . Therefore, L is independent of 3^' as well. We let 
1 h be the natural filtration of the process L t . 

The subordinate filtration 5Fj, the auxiliary subordinate bivariate process (X^,D^), and the 
default time Recall from our above discussion that & x , 5F D and 5F L correspond to the filtrations 
generated by the processes X, D, and L, respectively, and that the filtration J X,D = J^VJ 11 is the smallest 
filtration that contains ZF X and in which Q is an S^'^-stopping time. Similarly, we can define the filtration 
grG _ cfX,D y <jL ^ Q ^ g ^ e smallest filtration that contains J X,D and in which T t is an increasing family 
of 9 rG -stopping times. Then the subordinate filtration 5Fj is constructed by time-changing the filtration 5" G 
with the Levy subordinator T, i.e., Jt — 3j, . Observe that since T is an ^-stopping time, then 3^ is 
the filtration containing all of the information of the bivariate process (X, D) prior to the stopping time 
T t . Consequently, one can define the subordinate bivariate process (X^,D^) := {(X<r t , Dr t ),t > 0} by time- 
changing the bivariate process {X, D) with the subordinator T. This transformati on is call e d Bo chner's 
subordination due to work on subordination of semigroups an d their generators by iBochnerl (Il949h . The 
subordinate process (X^,D^) is a Feller J-semimartingale (see iMendoza- Arriaga and Linetskvl (|2012j)). 

From (|2.7p we observe that, before subordination, the indicator process D satisfies D t = 1 for all t > £• 
On the other hand, after subordination, the subordinate indicator process D* satisfies Df = 1 for all t > 
for which Tt > C (i- e -j Df = I/y t >£\). Therefore, in the credit-equity context one can define the default time 



by 



C cb = L C -, L a _=M{t:T t >a}. 



Certainly, is the first passage time process of T across the level £ and the identity {C^ < t} — {T t > C} 
holds (see section 5.1.1 in I Jeanblanc et aL (2009)). Hence, Df = I{r { >C} — I{c*<t} i s the default indicator 



process. The characterization of the subordinate process (X^,D^) is provided below in section [3l 

The stock price S and the scaling constant p. From Eq. (|2.1[) we observe that the dynamics of the 
stock price S can be described by means of the subordinate bivariate process (X^,D^). Indeed, the stock 
price S can be seen as a function f(t, Xf, Df) £ C([0, oo) x E) which is decomposed according to (|2.8|l with 
the payoff function u{t, x) — e pt x if no default occurs by time t > 0, and zero recovery v(t, x) = if the firm 
defaults prior to time t > 0. That is, 



S t = e pt Xf{l - Df). (2.11) 



The sc aling constant p is introduced to ensure that the asset S is an ^-martingale. As shown in lMendoza- Arriaga et al 



(l2010h . S will be a martingale if and only if p = <fi(—p) and p E 3, where the set 3 is defined in Eq. (|2.9p 



That is, assuming zero interest rates, K[St] < oo for every t > 0, and E[5t 2 l^tj = Sjj for every t\ <t%. From 
the previous condition in p we are free to choose any value of p as long as p G 3. Hence, from this point 
onward we assume that p, 6 3. Observe that the underlying assumption v(x) — is equivalent to modeling 
the stock price S under absolute priority, which means that the stock holder has zero recovery in the event 
of default. 



G 



3 Markovian and Semimartingale Characterization of the Subor- 
dinate Process (X*, D*) 

Before proceeding with the calculation of the quadratic variation of the log price process log S it is essential 
to describe the characteristics of the underlying stock process S in terms of the subordinate bivariate process 
(X^,D ^). For convenience, we summarize some of the key results of section 3 in lMendoza-Arriaga and Linetskv 
(2012 ) w ho give the Markovian and semimartin gale characterization of the process (X^,D^). We refer the 

reader to iMendoza-Arriaga and Linetskvl (|2012f) for the corresponding proofs. 

We begin by recalling some key results about subordination (in the sense oflBochnerl (I1949M of semigroups 
of operators in Banach spaces. The expression for the generator is due to IPhillipsI (1952). 

Theorem 3.1. (|Phillipsl f|l952h ~) Let T be a subordinator with Levy measure v, drift 7 > 0, Laplace 
exponent 4>(X), and transition function ■K t {ds). Let 7 be a strongly continuous contraction semigroup of 
linear operators on a Banach space 58 with infinitesimal generator A. 
(i) Define 

(3.1) 



7?f(x) = / 7 s f(x)n t (ds), t > 0, / £ B. 
J [0,00) 

Then 7^ — {7f,t > 0} is a strongly continuous contraction semigroup of linear operators on 58 called 
subordinate semigroup of V with respect to the subordinator T. 

(ii) Denote the infinitesimal generator of 7^ by A®. Then the domain of A is a core of A^ and 



A*f = jAf + 



(y s f - fMds), f £ Dom(A). 



(n, 



(Hi) Moreover, if 7 is a Feller semigroup on C([0, 00]), then the subordinate semigroup 7^ is also a Feller 
semigroup on C([0, 00]). 

Next, recall that if is not an absorbing boundary for the process X with diffusion coefficient a(x), drift 
b(x) and killing rate k(x), then the transition kernels P /3 (x, dy), j3 = 0, 1, of the semigroups 7@ have densities 
with respect to the Lebesgue measure, Pf(x,dy) = p^ (t,x,y)dy, where p@(t,x,y) are jointly continuous in 
t, x, y. This follows from the fact that any one-dim ensional d iffusio n h as a density with respec t to th e speed 
measure that is jointly continuous in t,x,y (cf. iMcKeanl ( 19561) or Borodin and Salminen ( 2002 ) p. 13). 
Und er our assumptions, the speed measure is absolutely continuous with respect to the Lebesgue measure 
(cf. iBorodin and Salminen ( 2002 ). p. 17) and, hence, the semigroups have densities with respect to the 



Lebesgue measure. For (0 = 0, the transition kernel is conservative, i.e., P®(x,E) — J E p°(t,x,y)dy = 1. 
For (3 = 1 the kernel is generally defective, i.e., Pf(x,E) = J E p 1 (t, x,y)dy < 1. While our diffusion is 
non- negative, for future convenience we extend the transition densities from E to K by setting p^(t, x, y) = 
for all y ^ E and for all x € E and t > 0. Then, the Markovian characterization of the subordinate bivariate 
process (X^,D^) can be obtained from Theorem 13. II as follows. 

Theorem 3.2. (Markovian characterization of (X^, D^)) (i) The bi-variate process (X^ , 
process with the Feller semigroup {7f,t > 0} acting on f £ C([0, 00] x {0, 1}) by: 



is a Feller 



7ff(x, d) = 7f\{x) + (1 - d)7p\u - v)(x), 



0,1/ 



where u(x) = f(x,0) £ C*([0,oo]), v(x) = f(x,l) £ C([0,oo]), and {7p v ,t > 0} and {7f ,t > 0} are 
Feller semigroups obtained by subordination in the sense of Bochner from Feller semigroups {7®,t > 0} and 
{7l,t>0}. 

(ii) The infinitesimal generator A^ of the Feller semigroup {7f,t > 0} has the following representation: 



A^f(x, d) = A*'°v(x) + (1 - d)A^\u - v)(x), u, v £ Dom(.A ), 



(3.2) 
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where A*^ , (5 € {0,1}, are generators of{yf' P ,t > 0}. 

(Hi) The generator A^'P has the following Levy-Khintchine-type representations with state- dependent 
coefficients: 



A+*f(x) = \ ia \x)f"{x) + b^(x)f'(x) - Pk+(x)f(x) 



+ 



(fix + y)- f{x) - I {lyl < iy yf'(x)y^(x, y)dy, f G Dom(^), 



where the state-dependent Levy density %^(x,y) is defined for all y ^ x by 



7T^(x,y)= p (s,x,x + y)v(d S ), (3.3) 

J(0,oo) 

and satisfies the integrability condition J m (\y\ 2 A l)ir^(x, y)dy < oo for each x G E (recall that we extended 
p@(t,x,y) iol by setting p@(t,x,y) = for y ^ E), the drift with respect to the truncation function h x4 ' (x) = 
xl {\x\<i} is given by, 

b^(x)=jb(x)+ [ (f yp^s,x,x+y)dy)u(ds), (3.4) 

J(0,oo) X J{\y\<l} ' 

and the killing rate is given by 

k*{x) = yk(x) + [ (l-pP(x,E))v(da), (3.5) 

J(0,oc) V ' 

where P@(x,E) = J E p l3 (s,x,y)dy. 

(iv) If fix, d) € Bom(A) (i.e., f is of the form (f278|) with u, v G Dom^ 1 )^ and (X^,D^) starts with 
Xq = x > and Dq = d G {0, 1}, then the process 



M{ := f(Xf, Df) - fix, d) - / A*f(X+, D*)ds 



is an 3 "-martingale. 



Now we turn our attention to the semimartingale characterization of the process (X*, D*) fsee lJacod and Shirvaevl 



(|2002f ). p. 76, for the definition of predictable characteristics of a semimartingale). 

Theorem 3.3. (Semimartingale characterization of (X^,D^)) (i) The bi-variate 3 '-semimartingale 
iX^, D®) has the following predictable characteristics. The predictable quadratic variation of the continuous 
local martingale component Xf' c is: 



Cf x * =1 f 'a 2 (Xf)ds 



(Cr D * = and C t x * D = since is purely discontinuous). The predictable process of finite variation 
associated with the truncation function (h x * (x) — xl{\ x \<iy, h ^ (x, d) = d) is: 

Bf = f b*> (Xf)ds, Bf = [\l - D*)k+(X})ds, (3.6) 

Jo Jo 

where the function b^'°(x) is defined in Eq. (|3.4p and k^(Xf) is defined in Eq. (|3.5[) . The compensator of the 

random measure [i x < D (u); dt, dy dz) — Ylu^-fAfx* d*)(u)/o/(«ax ( '(u) ^.{ds,dy dz) associated to the jumps 

of(X'^,D'f') is a predictable random measure on R + x (R 2 \{(0, 0)}): 

v x *' D \w;dt,dydz) = [ir*>°(x,y) - (I - d){^{x,y) - ^\x,y))] dyS Q (dz) 

+ (1 - d)jk{x)6oidy)6 1 (dz) + (1 - d)i^(x,y) - ^\x,y))dy 5 x {dz), 
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where n^'^(x,y) are the Levy densities defined in Eq. (I3.3[) with j3 = 0, 1, and S a is the Dirac measure 
charging a. 

[ii) The Levy-Ito canonical representation of X^ with respect to the truncation function h x (x) = 
xl{\ x \<i} is: 

X? = x + B?* +Xp c + f f yI {lv \< 1} L x \ds,dy)-v x \ds,dy)) + f f yl {lyl>1}f i x \d s ,dy), 

JO JR K ' JO JR 

where the compensator of the random measure fi x * (w; dt, dy) — ^ u I{ax*(w)^o}^(« ax*(u)) (^ s > dy) associated 
to the jumps of X^ is a predictable random measure on R + x (M\{0}): 

v x \uj;dt,dy) = Tr' t >' a (X?_,y)dydt. (3.7) 

(Hi) The Doob-Meyer decomposition of Df is: 

Df = Bf + Mf (3.8) 

with the martingale Mf = Df — and the predictable compensator B^ given in Eq. (|3.6[) . so that the 
^-intensity is Xf = (1 - Df)k^(Xf). 

Lastly, we formulate Ito formula for functions of the bi-variate process in the form convenient for our 
application. Since for p = <f>(— fi) the stock price process S is a martingale (hence, special) then it suf- 
fices to_jjresent L lJx^^ the case in which the process f(t,Xf,Df) is a special semimartingale 
(see. Jacod and Shiryaevl (|2002]), Definition 4.21, p.43). 

Theorem 3.4. (Ito Formula for (X*, D*)) Suppose (X*, D*) starts from xf = x > and D$ = d € 

{0, 1}. For any function f(t, x, d) — v(t, x) + (1 — d)(u(t, x) + v(t, x)) with u(t, x) and v{t, x) € C 1,2 (R+ x 
(0, oo)) (recall that zero is an unattainable boundary for the diffusion process X starting at x > 0), if 
f(t,Xf,Df) is a special semimartingale, ltd formula can be written in the following form: 



f(t,X?,Df) = f(0,x,d)+ / (d s +A*)f(s,X£,Df)ds + / d x f(s,Xt,Df)dXt' c 

Jo Jo 

+ f [ (v(s,Xt_+y)-v(s,Xt_j) {v x \xt-dsdy)-v x \xt-dsdy)). 

+ L L " ^ X °~ + V) ~ {U ~ ^ X ^ (1 " D s-^(Xt;dsdy) - u(xf_;d S dy)). 
\l-Di_){u-v){xt)dMt. 



where we introduced a random measure associated to those jumps of X^ that do not coincide with jumps of 

(i(U]d8,dy) =5Z I {A^( W )^0} I {A£.*M=0} 5 (t 1 ,AX*(c))( ds '^)' ( 3 - 9 ) 
u 

and its compensator measure 

v{u;ds,dy) = [^(Xf_,y) - (1 - Di_){^{xf_,y) - n^(Xt_,y))] dyds. (3.10) 

And where is the random measure associated to the jumps of X^ and v x4 ' is its compensator mea- 

sure (|3.7[) . The generator A® is given by Eq. (I3.2[) . 

Observe also that the continuous local martingale part that can be represented as xf' c = y^y J Q a(Xf)dW s 
where W is a Brownian motion, which will be convenient for our representation of the stock price process S. 
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4 Variance Swap Computation 



Due to the semimartingale characterization of the process (X^,D^) of section [3] we are now in position to 
provide the characterization of the stock price process S and to explicitly compute the value of variance swap 
rate K var = E([log5]t). First, let us recall from Eq. (I2.11[) that the stock price S can be seen as a function 
f(t, Xf, Df) 6 C([0, oo) x E) which is decomposed according to (|2.8[) with the payoff function u(t, x) = e pt x 
if no default occurs by time t > 0, and zero recovery v(t, x) = if the firm defaults prior to time t > 0. The 
following theorem formally characterizes the stock price process. 

Theorem 4.1. (Stock Price Process S) Let the stock price process S be specified in terms of the bivariate 
process (X* £>*) by the prescription S t = e pt (1 - Df)xf , where X$ = x > and D$ = d e {0,1}. 
Moreover, assume that the scaling factor p satisfies p = <f>{—p) where <f)(u) is the Laplace exponent (|2.10[) of 
the subordinator T , and where p £ J is the constant drift of the background process X (|2.2[) - (|2.3I) (the set 3 
is defined in Eq. (|2.9[) ). Then, the stock price process S is a martingale with canonical representation 



St — So 



V7 / a(S u )S u dW u + f y(l-Dt)(Kdsdy)-TT*' 1 (S s -,y)dyds)- f S u _dMjf, 
Jo Jo Jo 



where 7 > is the drift of the Levy subordinator T. The random measure \± corresponds to those jumps of 
X^ that do not coincide with jumps of the default indicator (see Eq. (|3.9p ). The Levy density ^^(x^y) 
is defined in Eq. (|3 . 3[) . W is a Brownian motion and is the martingale (|3.8p associated to . 



Proof. From the restrictions on p and p, the stock price St is a discounted martingale (see. lMendoza-Arriaga et al 
( 201dh . section 4). In the presence of a constant interest rate r > and dividend yield q > 0, the stock price 



can be decomposed as St — So + At + J S u —dM u , where Mt is a martingale and where At = / (r — q)S u -du 
is predictable, and hence, it is a special semimartingale. Then, the canonical representation of S follows 
from Ito's formula of Theorem 13.41 applied to the function f(t, x, d) defined by u(t, x) — e pt x and v(t, x) = 0. 
We further observe that the drift vanishes since 

(d s + A^)u(s, Xf, Df) = e ps Xf(l - Df) 

x 'p + 1 [p + k(Xf)] + (l/Xf) 



(0,00) 



(z - Xf)p x (w, Xf, z)dz)v{dw) 



jk(Xf ) - 


[ 




- f p 1 (w,Xf,z)d, 
Je 


e ps Xf(l - 


Df) 


P + 


J{0,oo) v 


e ps Xf(l - 


Df) 


P + 


1P+ 1 (e pw - 

J(0,oo) 


e ps Xf{l- 


Df)( P - 


^(-p)) = 0. 



((1/Xf) f zp 1 (w,Xf,z)dz-iy(dw) 



In the third equality we used the fact that E[e~ & k(x u )du X ^ = xe ^t ^ |Linetskvl (|2006h . Proposition 2.1). 
The fourth equality follows from the definition (|2.10[) of the Laplace exponent, which cancels from the 
condition p = <f>(—p). □ 

This canonical representation decomposes the stock price process S into a purely discontinuous martingale 
of jumps prior to default with the compensator measure (1 — D% _ )7r 1 ^ (xf_ , y)dydu (observe from Eq. (|3.10[) 
that (1 — D^_)D(du, dy) = (1 — Df l _)TT 1 '^(xf_, y)dydu), a continuous martingale component represented in 
terms of a Brownian motion, and a final jump to zero (the default term — J Q S u -dM!f). Clearly, the process 
S is a jump-diffusion process whenever 7 > 0, and a purely discontinuous process for 7 = 0. 

Next, we note that if the firm underlying S were to default at some time ^ in the interval [0,t], then 
the payoff of a traditional variance swap contract Ijl.ljl would be infinite. To account for this possibility, we 
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modify the floating leg of the VS so that it only accumulates quadratic variation prior to the default time 
Cj , That is, the long side of a VS, under our modified definition, has a payoff of 



[logS] tAC *_-K var = f (l-D£)d[logS] u -K va 
Jo 



(4.1) 



Notice that, for an asset that cannot default, our modified definition of a VS coincides with the traditional 
definition of a VS. Meanwhile, for an asset that can default, our modified definition of a VS is guaranteed 
to have a finite payoff, since the floating leg of the modified VS only accumulates quadratic variation up the 
time just prior to default. 

Using definition (|4.ip , the fair value of K var is the risk- neutral expectation of the floating leg 



K — 



(l-Dt)d[logS] u 



(4.2) 



An explicit expression for the right-hand side of (|4.2[) is given in the following theorem. 
Theorem 4.2. Let S be given by St = e pt xf(l — Df). Then the right-hand side of (|4.2[) is given by 



Ex [l {C *>u}7^W)]^+ / E. 



i {c * >u} J log 2 (l + -^-y^(xt , y)dy 



du. (4.3) 



Proof. In view of Eq. (I4.2[) , it suffices to first calculate the Levy-Ito canonical representation of the function 
f(t, x, d) = (1 — d) \og(e pt x), which corresponds to a zero-recovery log-contract on the stock price. That is, a 
contract that pays u(t, Xf) = \og(e pt Xf) = log S t if no default occurs by time t > 0, and zero otherwise (i.e., 
we set v(t, Xf) = 0). Hence, the canonical representation of the pre-default log-contract of S can be obtained 
by means of an application of ltd formula of Theorem 13.41 to the function f(t, x, d) = (1 — d) \og(e pt x), 

d(logS t ) =[(1 - Dt_)(dt+A^) \og(e pt xt)] dt + (l- Df_) ^fa(xf_)dW t 
+ (1 - Df_) f log (l + -^-) (Kdtdy) - 7T^(xt,y)dydt) 



(l-Dt)log(e pt xf_)dMf. 



(4.4) 



Observe that due to the default term (1 -Df_) log(e pt X f_)dAif = (l-Df_) \og(S t -)dMf the process jumps 
to zero at default, which is consistent with our selection of the function f(t, x, d) = (1 — d) \og(e pt x) that has 
zero-recovery in case of default. Consequently, it describes the pre-default dynamics of log St and prevents 
it from exploding at default time. From (|4.4|) it is straightforward to compute the differential d[logS]t of 
the pre-default dynamics of log St , 

d[\ogS] t = (1 - Df_) i^a 2 {Xf) + J^lo g 2 (l + Jl^y^{xt,y)dy^ dt 

+ (l - Df_) f log 2 (l + JU U(dtdy) - ^ 1 (xf_,y)dydt) 

+ (1 -Df_) log 2 (e pt xi)dDf. (4.5) 

Finally, multiplying flOJ by l-D$ = I {t < c * } , observing that (l-I>f_)(l-£)f) = (l-Df) and (l-Df)dDf = 
a.s., integrating over the interval [0, t], taking an expectation, and using the fact that the random measure 
(1 — Df_)[p,(dt dy) — n^' 1 (Xf_,y)dydt) is a martingale measure, one arrives at (14.31) . □ 



Next, we give an alternative formulation of the value of K var in terms of Feynman-Kac semigroups. 
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Proposition 4.3. Let S be given by S t = e pt Xf(l-Df) with X Q = x > and D% = d G {0, 1}. Also, let 7\ 
(resp., yf' 1 ) be the (resp., subordinate) Feynman-Kac semigroup defined in Eg. (|2.5I) (resp., Theorem ] 3. 2)) . 
Then, K var can be represented as follows, 

tf»«r = %*><>} 7 / {'?t 1 a 2 {x))du + 

I{C*>o} / ( / f{ s i ){x)v{ds)\du 1 (4.6) 

JO JO K J(0,oo) 7 

with 

f(s,y) = 3>J log 2 (y) - 21og(y)^ log(y) + log 2 (y)^l, y = X*. 

Proof. From Theorem l3.2l we observe that E[(l — Df)f(xf)] = (1 — d)y f' 1 f(x). Therefore, the first term of 
Eq. (|4.6[) follows immediately. From Proposition 32.5(m) in ISatol ( 19991 ). p. 215, we know that if ||CPi/(a;) — 
f(x)\\ = 0(t) as U 0, then / R f{y)n+(x, y)dy = f (0oo) (y s f(x)-f(x))v(ds). Since f(y) = \og 2 (y/x) = 0(\x- 
y\ 2 ) as y — > x, then to prove \\y\f(x) — /(a;)|| = 0(t) suffices to show that J E (y — x) 2 p 1 (t,x,y)dy = 0(t) as 
1 10. Indee d, the later hold s true since for an arbitrary e > 0, we have J E ^{\ x - y \<e}{y — x^p 1 (t, x, y)dy < Ct 
as t j (c£. iMcKearj (|l956jh Theorem 4.5). Therefore, 



iog 2 (i + ^ r y^(xt,y)dy 

u— 



log 2 (^-)p 1 (s,Xt,y)dy)u(ds) 
{y\ log 2 (X n _) - 21og(X u _)^ Iog(X„_) + log 2 (X u _)^l)^(d S ). 

(0,oo) 

The rest follows from observing that E[(l - D*)f(s, X%)\ = (1 - d)"?*' 1 f(s, x). □ 



5 Spectral Expansions 

In order for the results of sections[3]and|4]to be useful, we need a practical way to construct the FK semigroups 
{y\,t > 0} and {yp 1 ,t > 0} as well as the associated transition densities and p^' 1 ^, x,y) (recall 

that we had dropped the super-index a, since we only need the case of a = 1). Spectral theory, or more 
specifically, the theory of eigenfunction expansions, provides a straightforward method of constructing these 
operators and functions. Below, we review some useful results relating to eigenfunction expansions. A 
detailed description of the spectral Theorem for self-adjoint operators in a Hilbert space is given in appendix 

m 

Recall that the FK semigroup {y\,t > 0} has infinitesimal generator A 1 (ETo) . With A 1 we associate a 
scale density s and speed density m 

s(i) = ° p (- £ *) ■ m w = m ° p (£ 2 M dv ) ■ 

where the point xq is arbitrarily chosen in E = (0, oo). The generator A 1 (|2.6[) with domain 

Dom^ 1 ) = {/ € L 2 (E,m) : A 1 f G L 2 {E,m)}, 

is a self-adjoint operator in the Hilbert space IK = L 2 (E,m). H Therefore, we have spectral representations 
for the operators A 1 and g(A 1 ), where g is any Borel-measurable function. 

3 A 1 is dense in J{ implies that A 1 has a unique self-adjoint extension A 1 with Dom(yi 1 ) = J{. We will not distinguish 
between A 1 and its extension A 1 . 
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Let i{>\ and A be the generalized eigenfunctions/ values of — A 1 . Note that, since A 1 is the generator of a 
contraction semigroup {7\, t > 0}, the eigenvalues of — A 1 are non-negative. The operator "P\ can be written 
as y\ = e tA (for a general Banach space, when t he generator A 1 is unbounded the latter is understood as 
a strong limit via the Yosida approximation (see Pazv (1983), Corollary 3.5)). Thus, using (|B.2[) . for any 
/ € 5£ we have 



CA = (V>A,/) > 



W>A,/) 



^A(a;)/(a;)m(x)da;, (5.1) 



where ip\ indicates the complex conjugate of The notation X)a(' ' ' ) ^ s shorthand for 



£(• 



(■■■) du, 



where <Td(— A 1 ) and a c {— A 1 ) are the discrete and continuous portions of the spectrum of —A 1 respectively. 

Similarly, using the functional calculus of Theorem IB. H the subordinated semigroup CPf' defined in 
Eq. (|3.1[) can be obtained as, 



A 



ip\(x)f(x)m(x)dx, 



w here <t>(\) is the Lapla ce exponent of T, defined in Eq. (|2.10[) . One should mention that the recent book 
of ISchilling et~aD (l2010h is an excellent reference for Bochner subordination of semigroups (for example, the 
last result above is obtained from their Remark 12.4, p. 133). 



5.1 Uniformly convergence of the discrete spectrum 



In general, when the spectrum is discrete, the spectral expansion (|5.1[) of the semigroup T 1 (and hence, J* 1 ) 
leads to an infinite series. When the semigroup T 1 is of trace class, then it is possible to establish uniform 
convergence for the expansions as follows. Assume that the eigenvalues of —A 1 satisfy the condition 



£e~ At <oo, 



Vt > 0, 



(5.2) 



so that the FK s emigr oup 

t > 0} is trace class (see Section 7.2 of lDaviesI (|2007i) ). According to Theorem 
7.2.5 of lDaviesI (|2007h . if {7>l,t > 0} is trace class, then the eigenfunctions ip\(x) are continuous functions 
with the global estimate I^a^)! < e At / 2 y / p 1 (t, x, x)/m(x) for all t > 0. Setting f(x) = S y (x) in (|5.1|) yields 
the transition density p (i, x, y) of the FK semigroup. 



m(y)^2e 



-At 



■<Px(x)tp x (y). 



(5.3) 



The sum on the right-hand side of (15.31) converges uniformly in x and y on compacts. This ensures that, 
in addition to the L 2 convergence, the eigenfunction expansion (15.11) converges uniformly in x on compacts 
for all / € L 2 (E,m) and t > 0. This follows from the Cauchy-Schwarz bound for the expansion coefficients 
I c a| < y/ (/) /), the eigenfunction estimate, and the trace class condition (|5.2j) . 

In this case, the spectral expansion for the subordinated FK semigroup {J^ 1 , t > 0} can be obtained by 
conditioning on the subordinator T t . For any / € IK we have 

3>t' 1 f(x)=E ;c [?y(x)] =E[E X [9b t f{x)\T t ]] =]TE[ e ~ ATt ] caVa(x) = ^e-^CA^A^), (5.4) 

A A 
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where c\ = (ij)\,f) and cj>(\) is the Levy exponent of the subordinator T. If we assume that the Laplace 
exponent <f> is such that 

£y*< A >*<oo, VOO, 

A 

then the subordinated FK semigroup {Vf' 1 , t > 0} is trace class. If we further assume that the eigenfunctions 
tpx of the FK semigroup {?\ , t > 0} have a bound independent of A on each compact interval K = [a, b] C E 
(that is, if there exist constants Ck , which depend on the compact interval K but are independent of A, such 
that |-0a(x)| < Ck for all A £ adi—A 1 )) then, in addition to the L 2 convergence, the eigenfunction expansion 
of the subordinated FK semigroup (|5.4j) converges uniformly in x on compacts for all / G L 2 (E,m) and 
t > 0. As above, setting f(x) = S y (x) in (|5.4I) yields the transition density p^' 1 ^, 2, y) of the subordinated 
FK semigroup 

p^t, x, y) = m(y) e-^Hxix^xiy). (5.5) 

A 

The sum in (|5.5[) is uniformly convergent on compacts in x and y. Note that the semigroup corresponding 
to the JDCEV process of section [B~2l is of trace class. 



6 Examples 

In this section, we compute K var (|4.3p explicitly (up to an integral with respect to the Levy measure v of 
the subordinator T) , when the background Feller diffusion X is modeled as (i) a geometric Brownian motion 
with constant killing rate and (ii) a Jump-to-default Constant Elasticity of Variance process. 

6.1 Example: geometric Brownian motion with default 

Perhaps the most widely recognized non-negative diffusion in finance is the geometric Brownian motion 
process (GBM). Here, we consider GBM with a constant killing rate, which is a diffusion of the form (|2.2[1 - 
(I2.3|) with constant parameters k(x) = k > and cr(x) — a > (excuse the abuse of notation). The generator 
A 1 of the FK semigroup {V^ , t > 0} and the corresponding speed density m(x) are given by 

A 1 = \a 2 x 2 d 2 xx + (n + k)xd x -k, (6.1) 

m(z) = 4fexp(2£-logxV £ = - °-. (6.2) 

xa^ \ a a 2 



Most commonly, the FK transition density p 1 (i,x,?/) of GBM with default is written as 

_ e- kt ( -(\ogy-\ogx-{n + k-o 2 /2)tf \ 

P (t, x, y) — _ — . exp 

££rV27rf \ 2cr^t / 



(6-3) 



Due to the fact that A 1 (|6.1I) is a self-adjoint operator on the Hilbert space "K = L 2 (E,m), with m given 
by (|6.2p . the FK transition density p 1 (i,a;, y) also has an (generalized) eigenfunction expansion of the form 
(|5.3[) . The eigenfunctions and eigenvalues of — A 1 are given in the following Theorem. 

Theorem 6.1 (GBM Eigenvalues and Eigenfunctions). Let the operator A 1 be given by (|6.1[) . The spectrum 
of A 1 is purely continuous: o^-A 1 ) = a^A 1 ). The improper eigenfunctions of —A 1 and the corresponding 
improper eigenvalues are 

ipu(x) = ^/^exp ((iu -£)ilogaA , A w = i (u 2 + + k, cue (-00,00), (6.4) 

where £ is ywen m ()6.2j) . 
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Proof. One can check directly that the eigenfunctions and eigenvalues of (|6.4I) satisfy the improper eigenvalue 
equation — A tpu, — A^Vw an d the boundedness condition (IB. 31) . The orthogonality relation (ipu,ip^) = 
S(uj - v) follows by noting that ± J e-^ u ~ u ^ x dx = 8(u)-v). □ 

Remark 6.2. Transition density (16.31) can be obtained by writing eigenfunction expansion (|5.3p with the 
eigenfunctions and eigenvalues of (I6.4p . making a change of variables z = ^logy, and using the Fourier 
transform of a Gaussian density 

1 f -(z-b) 2 \ f 1 f. (h , 

■■ ex P Kit = / — — exp I i (b z)lo — du. 



V2vra 2 V 2fl2 / i(-oo,oo) 27r 

Remark 6.3. When the underlying S is given by (|2.1j) and the background diffusion X is modeled as a 
GBM with default, the at-the-money skew of the model-induced implied volatility surface is controlled by 
/i. For /i < jumps in S will be preferentially downward, causing a negative at-the-money skew. For \i > 0, 
jumps in S will be preferentially upward, causing a positive at-the-money skew. As the skew for equity 
options is typically negative, it makes sense to choose fi < 0. 

We are now in position to compute (I4.3[) when X is modeled as a GBM with default. 
Proposition 6.4. Let X be a GBM process with default as described above. Then we have (i) 

^ [l {c , >u}1 a 2 {Xt)\ du = 7 <J 2 ( ~ fc 6 j , 

and {ii) 

{c , >u} ^log 2 (l + ^y<f>'\Xt,y)dy du = a 2 ( ^^P^ j | 0oo e- ks (s + s 2 e)v(d S ). 



t 



'(o. 

Proof. See Appendix lA.il □ 
6.2 Example: Jump-to-Default constant elasticity of variance 

The Constant Elasticity of Variance (CEV) model of [Co3 (|l975h is a non-negative diffusion of the form 
([ ^ -([23 ]! . where = and 

ct(.t) = ax' 3 . 

Here, /3 < is the volatility elasticity parameter and a > is the volatility scale parameter. The specification 
f3 < is consistent with the leverage effect (volatility increases when the stock price falls). For f3 < the CEV 
process hits zero with positive probability. In particular, for /3 G [—1/2,0), the origin is an exit boundary. 
Fo r < —1/2 the orig i n is a regular boundary specified as a killing boundary. 

Carr and Linetskv (2006) extend the CEV model to include a possible jump-to-default. Their model is 



refereed to as jump-to-default CEV or, more succinctly, JDCEV. In the JDCEV framework, the jump to 
default has a killing rate which is an affine functions of the local variance 

k(x) =b + co- 2 (x) =b + ca 2 x 20 , 

where b > and c > 0. Although for all c > default may only occur through a jump from a positive value. 
When c > 1/2 the zero boundary is entrance for the JDCEV diffusion, and thus, the diffusion cannot reach 
zero from the interior of E. The majority of the expressions developed in this section hold for all c > 0. 
However, the credit-equity modeling framework developed in Sections\^^ works exclusively for the case in 
which c > 1/2 (i.e., the case in which zero is an entrance boundary). Therefore, one should keep in mind 
this restriction when applying the following more general results. 
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For a JDCEV diffusion, the generator A 1 of the FK semigroup CP 1 and the corresponding speed density 
are given by 

A 1 = \a 2 x 2fi + 2 8 2 xx + ([i + b + ca 2 x 2 ?) xd x -(b + ca 2 x 2fi ) . (6.5) 

m(x) = ^x 2c - 2 -^e X p(eAx-^), A = i^±M e = sign(/i + b). (6.6) 

a z a r\p\ 



The FK transition density p l {t,x,y) for the JDCEV diffusion was obtained bv lCarr and Linetskv ( 2006 ) 



. m(x)\fi + b\(xy)2 c e uu V 2 f A x~ 2 P +y~^ x \ r f A(xy)-P \ 
p 1 ^, x, y) = v ;|p \±Jt> exp -eA— y — Ait I v , \ , , 6.7) 

where l v is the modified Bessel function of order u, the constants A and e are given in (16.61) and 

2|P| [1Mb + < 

Due to the fact that the unique extension of A 1 is a self-adjoint operator in the Hilbert space = L 2 (E, m) 
with m(a;) given by (|6.6I) . the FK transition density (|6.7|) has an eigenfunction expansion of the form 
(15.31) . The eigenfunctions and eigen values of —A 1 are given in the following theorem, which is due to 



Mendoza-Arriaga and Linetskv ( 2010f h 



Theorem 6.5 (JDCEV Eigenvalues and Eigenfunctions). Let A 1 be given by (|6.5|) . When \fi + b\ ^ 0. 
the spectrum of A 1 is purely discrete: o~{A l ) — <Jd{A l ). The eigenfunctions of —A 1 and the corresponding 
eigenvalues are: 



Mx) = A»pJ in + W» + % exp + e)Ax- 2 ?) L^Ax^), (6-9) 

X n =oj(n-l)+X 1 , n = 1,2,3,- ■-, (6.10) 
where L v n are the generalized Laguerre polynomials, A and e are given in (|6.6p and v , u and \\ are given in 



Proof. One can verify directly that the eigenfunctions and eigenvalues (|6.9p - (|6.10p satisfy the (proper) eigen- 
values equation —A i/j n = X n ipn- Orthogonality of the eigenfunctio ns fiAn,, lAm) = fin.™, follows from the 
orthogonality relations of the generalized Laguerre polynomials (see lAbramowitz and Stegunl (|l972j ). pp. 
775) 



x a e- x Ll{x)L a m {x)dx = r( " + " + \ t , m . 

TV. 



□ 



Remark 6.6. The transition density (|6.7|) can be recovered from the eigenfunction expansion (|5.3I) with 



the ei genfunctions and eigenvalues defined in (|6.9I) - (|6.10I) by means of the Hillc-Hardy formula (see lErdelvi 



(119531) . p.189): 

t 1 ,Kia)Lm = exp (JJ±±M) j v ( ^ 

4^ r(n + v + 1) 11 w " w 1 - 1 V 1-i / 1 1-i 



n=0 

which is valid for all \t\ < 1, v > — 1 and a, b > 
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We are now equipped to compute K var (|4.3p when the background process X is a JDCEV diffusion. We shall 
focus specifically on the case fi + b < 0, since in this case, all of the relevant functions / are in L 2 (E,m) with 
m(x) given by (|6.6[) . Thus, we can compute all the necessary expectations explicitly using the eigenfunction 
expansion techniques of section [5] 

Before computing the expectations in (|4.3p it will be useful to give an analytical solution for the p-th 
moment of the stock price 



Ex 



(StY 



3 ppt E, 



3 ppt yM^p 



Proposition 6.7 (p-th Moment). Let the diffusion {X t ,t > 0} be a JDCEV process with parameters f3 < 0, 
a > 0, b > 0, and c > 0. Assume /u, + b < 0. Then, (i) for p > 2(f3 — c), the expected value of the function 
f(x) = x p is given by the eigenfunction expansion: 



E, 



(S t y] =e ppt Y^^ (Xn)t ~^n(x), 



where <fi(X) is the Laplace exponent of the subordinator T . The JDCEV eigenvalues X n and eigenf unctions 
ip n , are given in theorem ] 6. 51 and the expansion coefficients are given by: 



v p + 2c 



y/(n- iy.\fj, + b\T(v + n) 



p + 2c 



1,2,- 



(6.11) 



where (z) n = z(z— 1) • • • (z — n— 1) is the Pochhammer symbol. Also, (ii) the spectral expansion is uniformly 
convergent for allt > 0, absolutely convergent at t = 0, and uniformly convergent at t = if p > (/? + 1)/2 — c. 



Proof. The proof of part (i) is obtained from Lemma 3.1 and Proposition 3.1 in lMendoza-Arriaga and Linetskv 
(|2010l ). For part (ii) we note that the semigroup is of trace class. In addition, we note that for x G 
[a, b] C E, the eigenfunctions i() n (x) satisfy the bound I^An^)! < C/n 1 / 4 < C for s ome C < oo indepen- 
dent o f n although it may depend on the range [a, b] (see inequality (27a) on p. 54 of Nikiforov and Uvarovl 
()1988l )). Moreover, since the expansion coefficients satisfy the Cauchy- Schwartz bound, |c n | < v/ (/, /), 
then for any / £ L 2 ((0,oo),m) the spectral expansion of CPf' 1 / converges uniformly for all t > 0. That is, 
J2^Li e~ Xnt \c n ipn(x)\ < C\J(f, f) J2™=i e~ A "* converges uniformly for all t > 0. In Appendix IA.2I we show 
that the spectral expansion of y < ^' 1 x p also converges absolutely at t = 0. In addition, if p > (/? + l)/2 — c, 
then the spectral expansion converges uniformly at t = 0. □ 

Proposition 6.8. Let X be a JDCEV process with parameters /3 < 0, a > 0, b > 0, and c > 0. Assume 
/i + b < 0. Then we have 



I {C«>u}0' 2 (^) 



du 



ia 2 A^-m +1 r(c/\(3\)J2 



(l/(2|^|) + l) n _ 1 



-0(A„)t 



^ y/(n-l)\\» + b\r(v + n) 
Proof. The expectation in (|6.12|) can be written explicitly as 



0(A, 



(6.12) 



7a 



l { ^>u } (Xt) 20 



du 



7 a 2 ^-W+ 1 r(c/| / 3|) f ( 



(1/(2|/3|) + !)„_! 
y/(n-l)\\n + b\T(i/ + n) 



1p n (x))du, 
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where the last equality is due to Proposition [5I7J One should note that if 3/3 + 2c > 1 then the sum converges 
uniformly at t = and the integral can be done term by term. Otherwise, observe that: (a) the series inside 
the integral is absolutely convergent for alH > due to Proposition 16.71 (and continuous for all t > 0), and 
(b) the Laplace exponent <j>{\) is increasing. Then, we can conclude that the resulting series (|6.12[) is also 
absolutely convergent, and hence the exchange of sum and integral is justified (i.e., we integrate term by 
term with J -ds for some e > and then take the limit as e \, 0). □ 

Proposition 6.9. Let X be a JDCEV process with parameters (3 < 0, a > ; b > 0, and c > 0. Assume 
fi + b < and 2c — 2/3 > 1 . Then we have 



E, 



/ log 2 (i + ^-)^ 1 (xt,y)dy 



A 4 i<3 




OC OO 



EE' 

+ \{°> n=l m=l 



3 -0(A m )t 



du 



lp m (x) X 



(n-iy. ft2 

M + 6|r(^ + n) 




(m _ i)! (i - n) fe e,n-i(^ + fc + i)ei_ 1 (£gfi) 



6|r(z/ + m) 



(^)n-i r (^ + 

(n-1)! 



E 

fc=0 



r(i/ + jfc + i)fc! 




r(i/ + fc + i)fc! 



where 



(1 - a + v). 



^[{p-u + v) 



(l-a + u) n . 



in 



r(a)[(^(a)-£ 



p — a + ^ 



E 



^ (p - a + vf 



u(ds), (6.13) 

(6.14) 
(6.15) 



and ip(a) = T'(a)/T(a) (with no subscript on ip) is the Poly gamma function. 
Proof. We start by mentioning that the functions O* result from the integrals 



e'(«) 



z°-- x e- z log 5 (z)V n (z)dz, 



5 €{1,2}, 



Re(a) > 0, 



(6.16) 



which are available in lPrudnikov et al.l ()1983l) . Eq. 2.19.6.1 and 2.19.6.3, pp.469. The restriction 2c- 2^ > 1 
is imposed such that the sums converge absolutely at t — 0. Indeed, it is easy to show that for all x > we 
have | log(x)| < (1/x + x) and log 2 (a;) < (1/x + x). Moreover, since the expansion of (1/x + x) converges 
absolutely at t = due to Proposition 16. 7\ then each of the series also converge absolutely at t = 0. The 
rest of the proof consists of integrating term by term. Details are found in Appendix IA.3I □ 



7 Market Implied Levy Subordinators 

Note that the value of K var (14. 3|) depends on the drift 7 and Levy measure v of the Levy subordinator T. 
One could, of course, compute the value of K var by choosing a specific drift 7 and Levy measure v. However, 
this parametric approach would lead to a considerable amount of model misspecification risk, as there is no 
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guarantee that the chosen subordinator would induce European option prices consistent with those observed 
on the market. An alternative approach would be to use knowledge of (liquidly traded and efficiently priced) 
European call and put options to constrain one's choice of Levy subordinator. In this section we will show 
that, when the background diffusion X is fixed, the drift 7 and Levy measure v of the subordinator T can be 
obtained non-parametrically from the implied volatility smile of i-expiry European options. This approach 
greatly reduces model misspecification risk, as the obtained subordinator induces i-expiry option prices that 
are consistent with those observed on the market. 

Let S be described by (|2.1[) . We assume that background diffusion X is given, but that the drift 7 and 
Levy measure v of the Levy subordinator T are unknown. Denote by C(i, x; K) the price of a European 
call option with time to maturity t and strike price K. Note that the price of a call with strike price K can 
be obtained from the price of a put with the same strike through put-call parity. We assume the existence 
of Europ ean call options at all strikes K € (Q, 00). While calls at all strikes K G (0, 00) do not trade in 
practice. iBondarenkol ( 20031 ) shows how to estimate the value of call at any strike, given the value of calls at 
a discrete set of strikes. 

Let ps(t,x,y) be the transition density of S under the risk-neutral pricing measure 

p s (t, x, y)dy = P x [S t G dy] . 
Note that S t G dy if and only if (1 - Df )e pt xf G dy. Thus, 

Ps {t,x,y) = e-fpM&xrf) = e-^^r^A^d/), y' := ye""*. (7.1) 



As iBreeden and Litzenbergerl (|1978I ) show, the transition density ps(t,x,y) can be implied from a semi- 
infinite strip of call prices. We have 

C(t, x- K) = E x [(S t - K)+] = [ (y - K)+ Ps (t, x, y)dy. (7.2) 

JE 

Differentiating both sides of (|7.2[) twice with respect to K, and noting that d 2 iK (y — K) + = 5(y — K), one 
obtains 

d 2 KK C(t,x;K)= f d 2 KK (y-K) + p s (t,x,y)dy=p s (t,x,K). (7.3) 

JE 

Setting our the two expressions (|7.1|) and (|7.3|) for ps equal to each other yields 

d 2 KK C(t,x;K) = e-^m(K')J2e- HX)t Mx)^ x (K'), K' := K e - pt . (7.4) 

A 

Multiplying both sides of (|7.4j) by tpx' (K') an d integrating with respect to K, we obtain 

S& K C(t,x;K)$x(K')dK = e-#Y i e-+W t il>x{x) f ^ x {K')^ x ,{K')m(K')dK 

A Je 

= ^ e -« A )^W(^,^) = e-^'^W. (7.5) 

A 

Note that we have used dK' = er pt dK and {tpx,ipx') = S\,\>- We solve (I73|) for tf>(\) 

= -1 , f * KK c(t,^m*K) . (7 . 6) 

If there exists a Levy subordinator T independent of X, which is capable of generating the prices of call 
options on the market, then its Laplace exponent, evaluated at A is given by (17.61) . Thus, we refer to T 
with Laplace exponent (I7.6P as the market implied Levy subordinator. 



4 In a working paper, ICarr and Leel (,2006) obtain E \e ATt l using similar methods. The authors do not deal with Levy 
subordinators specifically. 
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Remark 7.1. It is worth noting that, although one can theoretically compute d] <K C(t, x; K) with call prices 
available at all strikes K € (0, oo), a more convenient expression for the integral in (|T. 6[) can be obtained by 
integrating by parts twice □ 



oo 

2 



d 2 KK C(t,x-K)^ x {K')dK = d K C{t 1 x-K)xj; x {K l ) q - C{t,x;K)d K M K ') 

C{t,x;K)d 2 KK ^ x {K')dK. 



We have the following limits 

lim d K C{t,x;K) = 0, lira C(i, x; K) = 0, 

K— >-oo K— >oo 

lim duCtt, x\ K) = — 1, lim C(i, a;: iO = x. 

A^O A-s-0 

Hence, we find 

/•OO 

5^C(f,i;if)i(]f')tf = ^(0) + i e-^ x Vx(0)+ / C(t,x; K)d 2 KK ^(K r )dK. (7.7) 

Jo 

Note that d x is a derivative with respect to the argument of 4>\ whereas d\ K is a derivative with respect to 
K. The advantage of using expression (|7.7p rather than the integral in (|7.6| is that the differential operators 
in (|7.7p act on the eigenfunction ipx rather than the call price C{t, x; K). Derivatives of ip\ can be computed 
analytically, whereas derivatives of call prices C(t, x; K) must be computed numerically from market data. 

Using (|7.6p one can obtain the value of 4>{X) for all A e a(— A 1 ). This information is sufficient for 
constructing the FK transition density p^' l (t, x, y) and the transition density ps(t,x,y) of S. However, to 
compute the value of K var we need the drift 7 and the Levy measure v of the subordinator T . As we show 
in the next two subsections, 7 and v(ds) can be obtained from limited knowledge of the map <f>. 

7.1 Obtaining 7 and v from 0(A): the compound Poisson case 

As noted in section [21 when the subordinator T is of the compound Poisson type, its Levy measure v can 
be written as the product of the net jump intensity a :— f((0, 00)) times the jump distribution F 

u{ds) = aF(ds). (7.8) 

In this scenario, the Levy-Kintchine formula (12.101) can be written 

0(A)= 7 A + aj^ {l-e- Xs )F(ds) = 7 A + a (l - F(A)) , (7.9) 
where we have defined F(X), the Laplace transform of the measure F(ds) 

poo 

F(\) := / e- Xs F(ds). 



The drift of the subordinator 7 and the net jump intensity a can now be obtained from </>(A) by taking the 
following limits 

d>(\) f f°° 1 — e~ Xs \ 
lim = lim 7 + a / F(ds) = 7, 



A— >oo V Jq A 

p 00 

lim (<f>(X) — 7 A) = lim a / (l - e _As ) F(ds) = a, 

A— »oc A— /n 



^Our thanks to Marco Avellaneda for pointing this out. 
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where we have used J °° F(ds) = 1. After obtaining 7 and a, we use (|7.9I) to solve for F 

F(A)=1 + 7A -^ (A) . (7.10) 
a 

Note, because er(— .A 1 ) C M+, equation (|7.10[) will not give us the value of F(X) for any A < (we know 
F(0) = 1). However, knowledge of F(X) for A < is not needed in order to uniquely determine F. To see 
this, we need the following theorem. 

Theorem 7.2 (Analyticity of the Laplace Transform). Let A € R. Let a(s) be a real, non-negative, non- 
decreasing function which satisfies a(0) =0 and is of bounded variation on [0, R] for every R > 0. If the 
integral 



2(A) := / e~ Xs a(ds), 
Jo 

converges for all A > X c , A c < 00 then 5(A) is analytic for all X > X c , and 

d n f°° 

—3(A) := j a (-sre-^a(ds). (7.11) 
Proof. See IWidderl (Il946l) . page 57, Theorem 5a. □ 



Remark 7.3. A function that is analytic in a domain D is uniquely determined over D by its values along 
a line segment in D. 

Remark 7.4. From (|7.11l) . it is clear that 3(A) is decreasing and convex. 

Let Dom(.F) = (A c ,oo) where A c G K. If the continuous spectrum of — A 1 is non-empty a c (— A 1 ) 7^ {0}, 
equation (|7.10|) gives us a map of F(X) on some interval / C Dom(F). The analytic extension of that 
map is unique and well-defined in throughout Dom(i ? ). If the the continuous spectrum of — A 1 is empty 
<T C (— A 1 ) = {0}, then (|7.10j) gives us a map of F(X) at a countably infinite number of points in Dom(F) (i.e., 
the proper eigenvalu es X„ of — A 1 ). In this case, the analytic extension of F(X) is still uniquely determined 
if F is Lipschitz (see Baumer and Neubrander |l994), Corollary 1.3) 



F(0) = and sup ^ (s)l < 00. 

t,s>0 \t — 8\ 

If F is not Lipschitz, the analytic continuation of F(X) is uniquely determined if we know the value of F(X) 
at equally spac ed intervals, i.e ., if, for n = 1, 2, 3, • • • , we know F(X n ) where X n = a + bn for some a > A c 



and b > fsee IWidderl (|1946l ). Theorem 6.2). Note that the eigenvalues (j6. 10[) of the JDCEV process are 
equally spaced. 

From a practical standpoint, one cannot evaluate (I7.10[) at an infinite number of A. Thus, in order to 
obtain F(ds) from (|7.10[) . one should seek to fit a positive, analytic, decreasing, convex function to a finite 
number of points of F(X). Upon doing this, one can use the inverse Laplace transform (Bromwich integral) 
to obtain F((0,s)) 



1 



(pC+ioo 



Here the constant C £ R is chosen so that the contour of integration lies to the right of all singularities 
of F (A). Another option for obtaining F(ds) is to numerically invert t he Laplace transform F(X ). For a 
survey of numerical techniques for Laplace inversion we refer the reader to lPavies and Martini ( 1979h and the 



references therein. Once F(ds) is obtained, the Levy measure v is given by (|7.8j) . In figure[TJ we graphically 
illustrate how to obtain 7, a and F(X) from knowledge of </>(A) at a discrete set of points. 
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7.2 Obtaining 7 and v from 0(A): the general case 



When the subordinates T is not of the compound Poisson type, its drift 7 can still be found using 

,00 j _ e _ As 

'0 



lim ^ = lim ( 7 

A— ► oo A A— >c 



A 



-v(ds) 



7- 



To obtain ^ we must introduce w(s) the toi/ 0/ f/ie Levy measure 



uj(s) :— v((s, 00)) = / ^(c£z) 

«/ (s,oo) 



Following iBertoinl (|2004f ). pp. 7, we note that 



0(A) 



7 = 



(0. 



-v{ds) 



e~ Xs w(s)ds =: w(A), 



(0,00) 



(7.12) 



where cD(A) is the Laplace transform of the function w(s). Depending on the nature of er(— .A 1 ), equation 
(|7.12j) either gives us a map of w(A) along a line segment / C Dom(cD), in which case the analytic extension 
of w(A) is unique, or (|7.12[) gives us the value of cD(A) at a countably infinite number of points. In this 
case, the analytic continuation of w(A) is uniquely determined if we know the value of w(A) at equally 
spaced intervals (the Lipschitz condition would not be satisfied for an infinite activity Levy process). From 
a practical standpoint, one may seek to fit an analytic, decreasing, convex function to a finite number of 
points of w(A). Upon doing this, one can obtain oj(s) from the Bromwich integral 



-1 />C+ioo 

oj(s) = / e As w(A)dA. 

2?n Jc-100 



Finally, one obtains v(ds) from 



v(ds) = —dio(s). 



8 Conclusion 

In this paper we model the price process of an und erlying S as a Feller diffusion ti me changed by a Levy sub- 
ordinator. This class of models, first developed in Mendoza-Arriaga et al. ( 2010t ). allows for the underlying 



to experience jumps with a state-dependent Levy measure, local stochastic volatility and a local stochastic 
default intensity. 

The contribution of this paper is two-fold. First, we show how to compute the price of a VS contract in the 
general Levy subordinated diffusion setting. Using our general formula, we perform specific VS computations 
when the background diffusion is modeled as (i) a GBM with default and (ii) a JDCEV process. Second, we 
show that the drift and Levy measure of the Levy subordinator that drives the price process can be obtained 
dirctly by observation of the i-expiry volatility smile. By using call and put prices to uniquely determine the 
Levy subordinator that drives the underlying price process, we reduce the risk of model misspecification. 

Our results open the door for future theoretical and empirical research. On the theoretical side, we plan 
to extend our results to price options whose payoff depends on both the realized variance and the terminal 
value of the underlying. On the empirical side, we plan to study the extraction of the Levy subordinator 
from market call and put data. 

Thanks 

The authors would like to thank Stephan Sturm and Vadim Linetsky for their helpful comments on this 
work. 
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A Proofs 



A.l Proof of Proposition 16.41 

Part (i) is a straightforward computation 

E» [h^>u}JV 2 (X£)] du = 1( j 2 f E x [I {c *> il} ] du = 7a 2 f e- ku du = 



1 -e 



-fci 



Part (m). Using Proposition 14.31 we have 



Hf*>«} 



Iog a (l + -S-V' 1 (^-.I/)d» 



du 



(0,oo) Jo 

2^ {^{x)\ogx) 



^og 2 {y)tp UJ (y)m(y)dy 



E\{x} 



lE\{x} 

^og(y)tp bJ (y)m{y)dy 



+ ?t' 1 (^(x)log 2 x) ■ U ^ u {y)m{y)dy ) } .LolurKl,). 



where 



yt l {^{x)\og n x) = (-ivYdZ (e-*^»M*j) . 
\og n (y)^Jy)m(y)dy = 2 ^a n ^ i n d%6(u) + 

IE 

Inserting (|A.2j) and (|A.3|) into (|A.lj) yields 

(S3) 



(0,oo) JO 

- 2 ( - 



(e-^^)"^(x)) • 2^a 2 - 1 / 2 * 2 d 2 ^ + if) 
+ (-io-) 2 9 2 (e-^^-^fic)) ' 2 V^~ 1/2 *(w + if) I duj du v{ds) 



2-1/2-2 



2 a 



(0,00) 



.'I - e -0(^)t 



1 - e-^-)* 



1 _ e -<A(A w )« 



S(ui + i^)dco v{ds). 



(A.l) 

(A.2) 
(A.3) 



(A.4) 



In the second equality we have integrated with respect to u and used integration by parts to move the 
off of the Dirac delta functions. Noting that 

if 9)" - Hfg'Y + 1 9" = f'g, iM(aO = VW^, 

we find 



A_i£ = k, 



= 2 v^ 2 - i/2 * 2 / k ( e - A -) • f 1 ^"^ m*) 

2 /l-e-tf*)* 



v(ds) 



4>{k) 



-ks 



(0,oo) 



;s (s + sV) u(ds). 
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A. 2 Computation of c n = (x p ,^ n ) from Proposition 16.71 

Part (ii). Consider the series U(x) = J2T=i u k( x ) f° r x E D C (0,oo). If for all x € D the function u^{x) 
satisfies the inequality |i tfc(x)| < dk. for fc = 1, 2 , • • • , where the series < oo, then the series {/(#) 

co nverges uniformly (see iPrudnikov et al. (1990). Section 1.3.4.3, p. 751). From the inequality (27a) on p. 54 
of Nikiforov and Uvarovl ( 19881) . we find that \ipi(x)\ < C/{n - l) 1 / 4 for some C < oo independent of n. 
Therefore, we have 



C 



\c n lpn(x)\ < 



(w„-l 



y/(n-l)\r(v + n)(n - l) 1 / 4 



To s how that the series J^,T d n converges, it is enough to show that for a large n we have log(d„)/ log(n) < — 1 
(see Prudnikov et al. ( 199dh Section 1.3.2.19, p. 751). Therefore, observe that 



d n 



n > 1. 



Thus 



lQg(^n) 

log(n) 



(„ _ l)^f+«/2-3/2 

(i/ + n)("+«)/2-i/4 ' 

) + n/2 - 3/2) log(ra - 1) - {(v + n)/2 - 1/4) log(i/ + n) 



i-p 

2|/3| 



< 



1-p 



n 
2 



log(n) 
v + n 1 
~~ 2 4 



log(n - 1) 
log(n) 



Note that log(n — l)/log(n) t 1- Moreover, it can be verified that the term inside the bracket is less than 
— 1 for all p > (0 + l)/2 — c (recall v = (2c + 1)/(2|/3|)). This shows uniform convergence at t = for 
p > (/3 + l)/2 — c. To show that the series is absolutely convergent at t = it suffices to show that 
lirrin-^oo dn+i /d n < 1 with d n = \c n ip n (x)\ (i.e., d'Alambert's test for convergence). Equivalently, J2n° 
converges if log(d„+i/d„) < as n — > 00. Hence, analyzing the asymptotic behavior and noticing that 
L n-i( z ) ~ e^z-^+^^in- 1)^/2-1/4 cos { 2 ^/( n _ i) z _ ^{2v + l)/4}/0F for n > 1, we find that 
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+ 7 +n- 



(iv + n)^")- 1 / 2 ' 



n > 1. 



Thus, we would like to test lim^-joo \og(d n +i/d n ) < 0, where n 3> 1. First observe that 



-(¥) = ( 



1 — p ^ 3\ , 
— ; — r H hn login ) — 

2 2 4/ 6V ' 
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v 

2 +n 



4 ) ], '" ! " 
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Then, making use of the approximation log(n + a) m log(n) + a/n — a 2 /(2n 2 ) + • • • , we find 



lim log ( 71+1 \ = lim 

n— >oc V d n ) n—toi 

which concludes the proof. 
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A. 3 Proof of Proposition 16.91 

First observe that from Proposition 14.31 we obtain, 

E, 



(0,oo) 

2 



l{C*>u} / \og 2 (l + -^ r y^ 1 (xt,y)dy du 

' Y, e ~ XnS {( \o g 2 ( y )MyMy)dy)[yt' 1 M^)] 

ly JE\{x} ' 



E\{x} 



]og(y)MvMv)dy) [^(log^J^Cx))] 



(y)m(y)dy) [O^ 1 ( log 2 (z) V>„ (x)) ] }d« i/(d. 



(A.5) 



where 



m—1 

oo „ 

^(log^C*)) = Yl e-^^di^mix), < m = ljog(y)My)^m(y)m(y)dy, 

m—1 

oo 

^(log 2 ^)^*)) = E e-^^dl^mix), d 2 n 

m—1 

Explicit representation of m and d^ iJTt are found using 
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Making the change of variable z = Ay~ 213 we obtain 

z"(log(z) - log(vl)) e - 2 i^ 1 ( i2: )L^_ 1 (z)^ 
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and 



1 / (n- l)!(m- 1)! 



4|/3| 2 V T{v + n)T(v + m) 



z v log 2 (z)e- z Ll_ 1 {z)L v m _ 1 {z)dz 



E 



log(A) / (n- l)!(m- 1)! 



2|/?| 2 y T{v + n)T{v + m) 



E 



z"\og{z)e--V n _ 1 (z)V m _ l (z)dz 



4|,«| 2 y r(z/ + n)r(z/ + m) 

Now, we use the identity f„ z u e~ z L l ^ l _ l {z)L^ n _ l {z)dz — 8 n>m T{v + n)/(n — 1)! (where 6„ jm is the Kro- 
necker delta) and the series expansion for the Generalized Laguerre polynomials (L^(z) = T(y + n + 
l)/n\ T,k=o(- n )kZ k /(T(is + k + l)fc!)) to obtain 



E 



and 

n-1 



- |y (n _ 1)!r(2 , + m) r(l/ + fc + m j e io & ^ e ^-1 w«* 



y (n - l)!r(i/ + m) ^ r(i/ + fc + ljfc! J E 



(n - l)T(u + to) ^ r(i/ + fc + l)fc! 



Similarly, we obtain 

\/(«- l)!|A« + &|r(i/ + n) 

where c„ is found by setting p = in (16.11| . Also, 
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l-2c 
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2G 



Substituting the above expressions into Eq. (|A.5|) and using the relation (16.161) for m , and 6^ m as well 
as equations (16. 14|) and (|6. 15|) . we arrive to the final expression (|6.13l) . 

B Spectral Theorem 

In this appendix we summarize the theory of self-adjoint operators acting on a H ilber t space . A d etailed 
exposition on this topic (including proofs) can be found in lReed and Simon d!98Ch andlRudin| (Il99lh . 



Let J{ be a Hilbert space with inner product (•,•)• A linear operator is a pair (Dom(.A), A) where Dom(.A) 
is a linear subset of J£ and A is a linear map A : Dom(.A) — » !K. The adjoint of an operator A is an operator 
A* such that (Af,g) = (f,A*g),Vf G Bom(A),g (E Dom(A*), where 

DomfTT) := {g G "K : 3 h G 5£ such that (A/, g) = (/, ft) V / G Dom(A)}. 

An operator (Dom(.A),.A) is said to be self-adjoint in !K if 

Dom(A) = Dom(A*), (AAflO = (/,-AflO V/, p G Dom(.A). 

Throughout this appendix, for any self-adjoint operator A . we wi l l assu me that Dom(.A) is a dense subset 



of ?£. A densely defined self-adjoint operator is closed f see iRudinl (|199lh . Theorem 13.9). 

Given a linear operator A, the resolvent set p(A) is defined as the set of A G C such that the mapping 
(A— Id A) is one-to-one and R\ := (A— Id A) -1 is continuous with Dom(i?A) = 5C- The operator i?> : IK — > J{ 
is called the resolvent. The spectrum o~(A) of an operator .A is defined as <r(.A) := C \ p(A.). We say that 
A G o~{A) is an eigenvalue of .A if there exists ^ G Dom(.A) such that the eigenvalue equation is satisfied 

yi^ = A^. (B.l) 

A function ip that solves (jB.ip is called an eigenfunction of A corresponding to A. The multiplicity of an 
eigenvalue A is the number of linearly independent eigenfunctions for which equation (|B.1|) is satisfied. The 
spectrum of an operator A can be decomposed into two disjoint sets called the discrete and essential E 
spectra: a (A) = crd(A) U a e (A). For a normal operator A, a number Ael belongs toaAA) if and only if A 
is an isolated point of cr(A) and A is an eigenvalue of finite multiplicity fsee iRudiri ( 199lh . Theorem 12.29). 



A projection-valued measure on the measure space (R, 23(R)) is a family of bounded linear operators 
{E(B),B G 23(R)} in IK that satisfies: 

1. E(Vt) = and E(R) = Id. 

2. £(£) is an orthogonal projection. That is, E 2 (B) = E(B) and E(B) is self-adjoint: E*(B) = E(B). 

3. E(AnB) = E(A)E{B). 

4. If B = |J^ 1 B, and B, n Bj ■ = for i ^ j then £J(B) = lim^oo £)* =1 J5(S 3 -) S where the limit is in the 
strong operator topology. 

5. For every in /, g G K the set function /i/ i£) (B) := (f,E(B)g) is a complex measure on 23(R). 

Theorem B.l (Spectral Representation Theorem). There is a one-to-one correspondence between self- 
adjoint operators A and projection-valued measures E on 3~C, the correspondence being given by 

A = f XE(dX). 

If <?(•) is a B orel function on R then 

9(A) = f 9(X)E(dX), Dom(ff(.A)) = {/ G IK : / | 5 (A)|V/,/(dA) < oo}. (B.2) 

6 The essential spectrum may be further decomposed into the continuous spec trum and the residual spectrum. It can be 
shown that the residual spectrum of an ordinary differential operator is empty (see [Roach (1982), page 184). 
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Proof. See lRudinl (|l99lh Theorems 12.21 and 13.33. □ 

As a practical matter, if A is a differential operator acting on a Hilbert space L 2 (I, m(x)dx), where / is an 
interval with endpoints I < r, then the operators defined by (|B.2[) can be constructed by solving the proper 
and improper eigenvalue problems 

proper: Atp n = A„ ip n , A„ G (Td(A), ip n G JC, 

improper: .A^ w = A w t/^, A w G cr e (A), tp^ </ Oi. 

For the improper eigenvalue problem one extends the domain of A to include functions all functions / for 
which the following boundedness conditions are satisfied 

lim |/(a:)| 2 m(x) < oo, lim |/(a;)| 2 m(a;) < oo, (B.3) 

We will use Latin subscripts (e.g., I, m, n) to denote proper eigenfunctions/ values and Greek subscripts (e.g., 
w, fi) to denote improper eigenfunctions/ values. When we do not wish to distinguish between proper and 
improper eigenfunctions/ values we will write ip\ and A with no subscript. We refer to ip\ and A as generalized 
eigenfunctions / values . 

After normalizing, the proper and improper eigenfunctions A satisfy the following orthogonality relations 

(ipn, 4>m) = 8 n ,m, (lpu>, Vv) = <K W _ M)> (lpn,1pu) = 0- 

The operator g(A) in (|B.2I) is constructed as follows (see lHanson and Yakovlevl (|2002[ ). section 5.3.2) 

g(A)f = Y f 9W{i>xJ)^x 
x 

: = X! 9(K) (i>n, f) i>n + / g(^u)(lpu,f)lpujduj- 
X n ea d (A) J\„e<re(A) 

It is not always easy to evaluate divergent integrals of the form {ipx^ipy) and verify that they are in fact 
delta functions 8{\ ~ X') . A method for d irectly obtaining properly normalised improper eigenfunctions can 
be found on pag e 238 of lFriedmanl (|l956l ). 



7 The term "improper" is used because the improper eigenvalues A ^ cr d (yi) and the improper eigenfunctions Tp\ ^ Jf since 

(i>u,ipu) = oo. 
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Figure 1: A graphical illustration of how to obtain the drift 7, the jump intensity a and the Laplace transform 
F(X) of the jump distribution F(ds) for a subordinator of the compound Poisson type from a discrete set of 
values of the subordinator's Laplace exponent <^>(A). Top left: a plot of 4>(X) for a Levy subordinator with 
exponentially distributed jumps v(ds) = ar)e~ r ' s ds. Top right: a plot of </>(A)/A for the same subordinator. 
The level of the solid line corresponds to drift of the subordinator 7. Bottom Left: a plot of </>(A) — 7A for 
the same subordinator. The level of the solid line corresponds to the net jump intensity a. Bottom right: 
a plot of F(X) and its analytic continuation for values of A < 0. In all four plots we use parameter values 
7 = 1/100, a = 1 and 77 = 10. 
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